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a Yukawa correction with a strength of 10“^ times gravity and length scale of 100 AU is within 
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gravity including many theories motivated by dark matter or dark energy. 
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In recent times a great deal of attention has been de¬ 
voted to searching for modifications of Newtonian gravity 
at or below mm distances [T]. This possibility is moti¬ 
vated for example by theories with large extra dimen¬ 
sions or supersymmetric theories with light moduli mg. 
At the opposite extreme of the distance spectrum, devia¬ 
tions from Newtonian gravity at very long distance scales 
may be related to astrophysical and cosmological prob¬ 
lems in modern physics |g. The cosmological constant 
problem is one of the deepest and most enduring myster¬ 
ies in modern physics. There have been many theories 
which attempt to modify gravity on long distance scales, 
for example Dvali-Gabadadze-Porrati (DGP) gravity or 
theories of massive gravity jZHig. In these theories large 
effects are predicted to arise only on cosmological length 
scales. However, smaller effects may be seen at shorter, 
experimentally accessible length scales niiiT]. Eurther, 
theories that replace dark matter with modifications of 
Newtonian gravity, such as modified Newtonian dynam¬ 
ics (MOND) [iHl [TO] may also lead to observable effects 
on long, yet still accessible scales [20]. The enduring neg¬ 
ative results [TOHTO] in the search of dark matter in the 
form of weakly interacting massive particles (WIMPs) 
may indeed suggest a broader approach to this important 
problem. Although theories with long-distance modifica¬ 
tions of gravity have historically been troubled by ghosts, 
discontinuities, and other theoretical difficulties, recent 
progress has been made in alleviating these concerns; 
for a review, see Ref. [22] ■ For these reasons, we focus 
on the directly observable consequences of long-distance 
modifications of gravity and estimate the potential of a 
dedicated space mission to measure modifications of the 
gravitational force of the Sun out to 100 AU. This is the 
largest distance accessible for direct measurements in a 
practical amount of time using known technology and 
would provide a powerful test of long-distance modifica¬ 
tions of gravity. 
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One way to parametrize possible deviations from the 
1 behavior of gravity is by introducing a new Yukawa 
force with charge proportional to mass, so that the effec¬ 
tive gravitational potential can be written as 


^{R) = 


GM 

nr 


1 + ae 


-R/A 


( 1 ) 


where G is Newton’s constant, M is the source mass, and 
R is the distance from the source. The new Yukawa in¬ 
teraction then has a strength a relative to gravity and a 
characteristic length scale A. In this framework, experi¬ 
ments measure (or constrain) the dimensionless parame¬ 
ter a as a function of the distance scale A. A summary of 
current limits on the magnitude of a is shown in Fig. 

Qualitatively, the sensitivity of experimental measure¬ 
ments to new Yukawa forces of strength a improves as 
A becomes substantially larger than Earth-scale inhomo¬ 
geneities and reaches a level < 10“^° for A ^ i?Moon 
(the radius of the lunar orbit) owing to laser lunar rang¬ 
ing measurements [6]. Beyond such distance scales, the 
orbital mechanics of the planets in the Solar System 
forms the most stringent test of new Yukawa contribu¬ 
tions to gravity [2^. The Pioneer 10 and 11 spacecraft 
provided an alternative way to measure a at very long 
scales A ^ 10-100 AU as they receded from the Sun over 
a period of over 30 years. Indeed for some time the anal¬ 
ysis of the Pioneer data appeared to indicate an anoma¬ 
lous acceleration towards the Sun of 10“^° m/s 2 Enum, 
resulting in |q;| ~ 1.7 x 10“^ at A ~ 10^^ m. A sub¬ 
sequent, more careful analysis, however, attributed this 
effect to instrumental systematics |28]. The many un¬ 
certainties resulting from spacecraft designed for other 
purposes make the upper limit on |a| not competitive 
with that computed from planetary dynamics. 

As will be discussed (and as shown in Fig. [^, a spe¬ 
cially optimized space mission offers the opportunity 
of improving sensitivity to a new Solar System-scale 
Yukawa force by at least two orders of magnitude with re¬ 
spect to the present state of the art. We show that such a 
mission can be relatively simple and would use only well- 
understood and tested technology. In addition, measure¬ 
ments directly done along a trajectory actually reaching 
100 AU constitute a test that is fully model independent 
and would detect anomalies that are not well-described 
by a Yukawa term. This is in contrast with measure¬ 
ments using planetary motion that rely on the specific 
Yukawa parametrization to extract deviations at scales 
A different from the orbital radius. While the Yukawa 
potential is quite generic, there are interesting modifica¬ 
tions of gravity that produce truly long-distance effects 
that are not accurately described by a Yukawa potential. 
Reliance on the Yukawa functional form is a limitation 
of the previous searches in Ref. [5S]. In a strict sense, 
no data exists beyond the distance of Jupiter except for 
the limits set by the Pioneer spacecraft. The mission 
proposed here would improve these direct limits by over 
four orders of magnitude by carefully controlling system¬ 
atic effects and performing direct measurements along 


A(m) 



FIG. 1. Existing 2a experimental limits on new Yukawa forces 
with strength |a| relative to gravity as function of the scale 
A. The grey region, adapted from [1], is the current state of 
the art. The dashed line is the size of the “Pioneer anomaly,” 
which can be interpreted as a limit set by the trajectory of 
the Pioneer spacecraft. The dotted curve corresponds to the 
expected sensitivity of the experiment proposed here. 


the journey to 100 AU. 

II. EXPERIMENTAL CONCEPT 

Two guiding principles are key in the design of a 
space mission optimized for a sensitive search for devi¬ 
ations from the 1/R^ law at large distances: the mis¬ 
sion should have a reasonable duration and the space¬ 
craft should be designed in such a way as to minimize 
the non-gravitational interactions on the body whose ac¬ 
celeration is being measured. Both issues were far from 
ideal in the Pioneer missions that were designed for the 
exploration of the outer Solar System. Pioneer 10 took 
37 years to reach 100 AU and substantial systematic un¬ 
certainties on the measurement of the acceleration vector 
a(i?) occurred due to thruster leakages, drag produced by 
interplanetary dust and solar wind, and recoils against 
various forms of radiation emitted by the spacecraft. 

The concept discussed here is based on a low mass 
(M = 200 kg) spacecraft propelled by a heavy rocket. Af¬ 
ter a series of planetary flybys designed to gain speed, the 
spacecraft would then coast, while performing the mea¬ 
surements. For concreteness, we assume that the coast 
phase takes the spacecraft from ^ 1-100 AU, consistent 
with a series of flybys via Jupiter, Mars, and finally. 
Earth. The simple model here conservatively assumes 
that the position and velocity of the spacecraft is mea¬ 
sured once every ~ 3 weeks. Likely, Deep Space Network 
(DSN) ranging will be available for measurement more 
often than this assumption. A preliminary and generic 
flight time calculation is consistent with a total coast 
time to 100 AU of seven years [H]. The sensitivity of 
the experiment depends somewhat on the details of the 
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DF feedback thrusters 



FIG. 2. Conceptual view of the drag-free (DF) and relay 
spacecraft. The figure is not to scale, and the DF craft is 
shown in cross-section. The DF craft uses micro-thrusters to 
remain centered on the free-falling PM while rotating in a 
plane perpendicular to the Sun to reduce thermal and grav¬ 
itational backgrounds. The relay craft carries a high-gain 
antenna, always pointing towards Earth, for DSN ranging 
and communication. Communication and ranging between 
the two spacecraft are performed by small omnidirectional 
transponders. 


coasting trajectory; for most of this paper we consider 
a polar trajectory, perpendicular to the ecliptic plane of 
the Solar System. This trajectory reduces the impact 
of the Kuiper Belt’s highly uncertain gravitational pull. 
Other trajectories are possible, but have different sys- 
tematics to consider; we discuss the choice of spacecraft 
trajectories in Sec. |IIB| 


A. Spacecraft 

Central to the spacecraft design is the use of the drag- 
free (DF) technique [3^, whereby a feedback system 
“flies” the spacecraft around a “proof mass” (PM) that, 
to within a very high degree of accuracy, is subject only 
to gravity. The PM is stowed and inactive during the 
initial maneuvering and flyby phases and is only released 
and tracked during the coast. The DF technique was de¬ 
veloped in the ’60s and initially tested on the U.S. Navy 
TRIAD spacecraft m- More recently, Gravity Probe 
B [32] used the DF technique in a challenging configura¬ 
tion, where each of the quartz rotors of the gyroscopes 
at the heart of the experiment were also used as PMs. 
In 2009, ESA launched the Gravity Field and Steady- 
State Ocean Girculation Explorer (GOCE) [33] that uti¬ 
lized the DF design to map Earth’s gravitational held. 
The Laser Interferometer Space Antenna (LISA) [SI] also 
plans to use the DF technique to establish a highly ac¬ 
curate geodesic network between PMs located on differ¬ 
ent (and widely separated) spacecraft. As will become 
clear, the DF system discussed here is conceptually sim¬ 
pler than those of Gravity Probe B and LISA, although 
new and unique challenges are presented by the long du¬ 
ration of the hight and the needs of telemetry and ranging 
over very large distances. 


Parameter 

Value 

DF Spacecraft Mass (M) 

200 kg 

DF Spacecraft Radius (rs) 

1 m 

Experiment Duration (t) 

7 yr 

PM Mass (to) 

10 kg 

PM Radius 

5 cm 

DF Cavity Radius (rq) 

10 cm 

Required Thruster Bandwidth 

10-2 Hz 

Optical Sensing Deadband (d) 

10 /im 

Minimum Correction Period 

100 s 

Time Between DSN Measurements 

2 X 10® s 

PM Discharge Period 

2 X 10® s 

Microthruster Fuel Mass 

< 50 g 

Angular Velocity {uj/2tt) 

0.1 Hz 

Spacecraft Velocity, Radial {v) 

14 AU/yr (initial) 

Trailing Spacecraft Distance 

10 km 

RTG Thermal Power Requirement 

1 kW 


TABLE I. Summary of relevant parameters, as outlined in 
the text. The PM mass assumes the use of platinum for its 
construction, and the fuel mass required includes only the 
amount necessary for operating the DF system. The RTG 
thermal power refers to the trailing spacecraft and would re¬ 
sult in ~ 50 W of electrical power. The drag-free craft may 
use a substantially smaller RTG or, maybe more likely, re¬ 
ceive power transmitted optically from the relay craft. The 
symbols in parentheses match those used in Table [Hj 


In our model we assume that measurements of the 
1/R^ law can be carried out over about two orders of 
magnitude in distance, from ~ 1 AU to 100 AU. While 
conventional chemical engines would be employed for the 
flyby phase, such engines would be jettisoned at the 
beginning of the coast when micro-thrusters optimized 
to assist the DF system over the long duration of the 
flight would take over. Standard telemetry using NASA’s 
Deep Space Network (DSN) [35] would provide range and 
Doppler data at well known times, from which a(i?) can 
be derived. 

A conceptual sketch of the spacecraft in the coasting 
configuration is shown in Fig. The ultimate perfor¬ 
mance of a DF system is limited by the interactions be¬ 
tween the spacecraft (shell) and the PM. These include 
thermal gradients on the inner surface of the shell, resid¬ 
ual gas from spacecraft components, and gravitational 
forces from asymmetries in the shell. The thermal gra¬ 
dients can be reduced with appropriate insulation, and 
the residual gas can be reduced by appropriate choices of 
materials. However, gravity cannot be shielded, so great 
care has to be taken to build the spacecraft symmetrically 
around the cavity hosting the PM. In order to minimize 
the residual gravitational interaction of the spacecraft on 
the PM at the maximum offset allowed by the feedback 
system, it is advantageous to design the gap between the 
PM and the inside of the cavity to be larger than found in 
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existing DF implementations. In the conceptual design 
discussed here, the cavity and PM are spherical, with 
the cavity radius rc = 10 cm or more and PM radius 
I"PM = 5 cm. These dimensions are not the result of a 
careful optimization but are derived in analogy with pre¬ 
vious systems [Ml Elj with the constraint of a modest- 
size spacecraft. With such a large gap it is likely that 
optical ranging would be appropriate for measuring the 
position and motion of the PM in its housing. A (prob¬ 
ably different) optical system would also be a candidate 
for initializing the DF system after the PM is released 
at the beginning of the coast. Ultraviolet light-emitting 
diodes (LEDs) inside the cavity would allow for periodic 
discharging of the PM, which accumulates charge over 
time due to interactions with cosmic rays, every few days 

[Miisn]- 

The choice of thrusters to be used for the DF sys¬ 
tem is constrained by considerations of reliability, fuel 
endurance and bandwidth. In addition, care must be 
taken to ensure that fuel consumption does not change 
the position of the mass center of the spacecraft, produc¬ 
ing an anomalous acceleration of the PM. The bandwidth 
/ and total impulse J requirements can be estimated by 
imposing a deadband of 10 /rm on the centering of the PM 
with respect to the spacecraft center of mass. This results 
in / ~ 10“^ Hz and J = 2 x 10“^ Ns. Electromagnetic 
thrusters provide high specific impulse, proportional con¬ 
trol and, in the case of field emission electric propulsion 
thrusters (FEEPs) [40], bandwidths in excess of 1 Hz. 
FEEPs emit and electromagnetically accelerate Cs, Rb, 
or In ions, providing thrusts of 10“^-10“^ N with the req¬ 
uisite bandwidth m- Using the specifications in Ref. [M] 
for a cesium-based FEEP, the calculated fuel consump¬ 
tion is < 50 g for the duration of the mission. However, 
additional fuel is likely needed for attitude control and 
for adjustments to the spacecraft rotation. While it is 
not clear if FEEPs with the required reliability will be 
available, alternative schemes with nominally similar per¬ 
formance include laser ablative microthrusters m and 
other types of electromagnetic microthrusters, currently 
under development |42] . 

In order to further reduce the effect of interactions with 
the spacecraft on the PM, rotation of the shell around 
the PM can be imposed during the coast phase. Such 
a rotation spectrally shifts all disturbances produced in¬ 
ternally to the spacecraft in the plane perpendicular to 
the rotation axis. Since, to first order, only accelerations 
along the coasting direction are important, the plane of 
rotation is chosen to contain the Sun. For the purposes 
of preliminary calculations, we assume a rotation rate of 
^0.1 Hz, which is sufficient to reduce many nuisance 
interactions to negligible levels while remaining easily 
achievable with onboard thrusters. 

Due to the rotation of the probe and to the need for 
a high gain directional antenna for communications from 
deep space, we envisage the mission to include a second, 
trailing spacecraft approximately 10 km behind the DF 
probe, to relay telemetry to and from Earth. The 10 km 


distance renders the gravitational coupling between the 
trailing spacecraft and the PM negligible, while being 
short enough to allow for simple omnidirectional com¬ 
munications. The DF and relay crafts would be docked 
together during the maneuvering phase of the mission, 
and then separate and take their relative positions at the 
beginning of the coast. The relay spacecraft does not 
need to accurately maintain its position relative to the 
DF probe, as long as their relative position and velocity is 
constantly measured with sufficient precision. In the rest 
of this paper we will assume that all such ranging and 
velocity measurements are the combination of the two 
segments: Earth-relay craft and relay craft-DF craft. 

Radioisotope thermal generator (s) (RTGs) will have 
to be used to provide power to the spacecraft because 
of the large heliocentric distance reached. We expect 
a conventional RTG rated for I kW thermal and SOW 
electrical power to be sufficient for the relay craft. A 
smaller RTG or a system to optically transmit power can 
be used for the drag free probe. Optical transmission of 
power may benefit from a flight formation in which the 
two crafts are separated in the plane orthogonal to the 
Sun-PM axis, so that the center of rotation of the drag- 
free probe is visible from the relay craft. 

The relevant parameters of the proposed spacecraft are 
summarized in Tabled 


B. Trajectory 

The spacecraft trajectory is a key factor in determining 
the sensitivity of the experiment. The trajectory deter¬ 
mines the effects of any new physics as well as Solar Sys¬ 
tem backgrounds. For instance, trajectories that remain 
quite close to the Sun (say, in a bound, nearly circular 
orbit) have larger effects from new physics, but cannot 
discriminate between new physics and backgrounds—in 
this limit, all effects are, to first order, perihelion preces¬ 
sions. Such experiments can set excellent limits with null 
results (lunar laser ranging (LLR), for example, sets the 
strongest limit on |a| of any existing experiment) but are 
unable to directly confirm any signal they may observe. 
For that reason, we focus on trajectories that traverse a 
wide range of distance scales, making possible the iden¬ 
tification of a signal based on its radial dependence. 

The Pioneer and Voyager spacecraft followed such tra¬ 
jectories, close to the ecliptic plane of the Solar System 
in order to conduct several planetary flybys and observa¬ 
tions. However, a trajectory in the ecliptic plane passes 
through the Kuiper Belt, which has a very poorly con¬ 
strained mass distribution. The systematic uncertainty 
introduced by the Kuiper Belt in this case is found to be 
substantial. A polar trajectory, coasting perpendicular to 
the ecliptic plane, reduces the effect of the Kuiper Belt 
while also allowing its mass distribution to be ht with a 
small number of parameters (see Appendix . We note 
here that future experimental constraints on the mass 
distribution of the Kuiper Belt may turn the ecliptic tra- 
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jectory into a competitive option, but that possibility is 
not considered further in this work. 

Reference [53] shows that spacecraft velocities of up 
to 14 AU/yr are achievable with a sequence of planetary 
flybys. Achieving these high speeds requires careful tra¬ 
jectory designs which we do not attempt here. Instead 
we present an estimate of sensitivity using a representa¬ 
tive trajectory whose final flyby is around Earth, out of 
the ecliptic, and at an initial speed of 14 AU/yr, reaching 
77 AU from the Sun in seven years. 


III. ANALYSIS OF SYSTEMATICS 

The fundamental limit to the sensitivity of the experi¬ 
ment is the accuracy with which the position of the PM 
is measured via DSN ranging. However, there are nu¬ 
merous other sources of error, resulting from external 
forces which act on the PM. Although the drag-free sys¬ 
tem eliminates the largest external forces on the PM, 
forces which act directly on the PM (either external or 
caused by the shell itself) will still affect the net motion 
of the spacecraft. These forces must be minimized and 
well-constrained in order to be able to positively identify 
a signal. 

Some forces (especially those from the shell itself) can 
be made sufficiently small so as to be negligible compared 
to the ranging measurement uncertainty. Requiring that 
the force lead to a displacement of less than 1 m after 
the full 7 yr coast of the spacecraft translates into a con¬ 
straint on the design of the shell and the precision of its 
construction. External forces, such as the gravitational 
forces from objects in the Solar System, cannot be engi¬ 
neered away. Instead they must be precisely modeled and 
subtracted from the motion of the PM. We find that the 
masses and positions of the planets are sufficiently well- 
known to model and subtract their effect [331 031 • The 
Kuiper Belt, despite its small total mass, has a mostly 
unknown mass distribution and cannot be subtracted or 
fit by a single parameter; we find that for a polar tra¬ 
jectory, three fit parameters are required, corresponding 
to the first three terms in the expansion of the Kuiper 
Belt’s gravitational potential discussed in Appendix [A] 

The sources of error and resulting design constraints 
are summarized in Table jlT] and fully detailed below. 

The forces acting on the PM that we wish to min¬ 
imize fall into a few categories based on the resulting 
displacement of the PM. We find that constant, R~^^ 
rotating, and random forces all come into play. A con¬ 
stant force leads to a displacement \aT'^ after a coast 
of time T. A force which falls off with distance (e.g. as 
i?“^) is most important at the beginning of the coast, 
where it builds up some additional velocity, which then 
adds to the displacement only linearly; the displacement 
at the end of the coast is roughly aoRoT/v, where oq 
is the acceleration at the starting position i?o — 1 AU 
of the coast, and u ~ 11 AU/yr is the average velocity 
of the spacecraft. Since the shell will be rotating with 


angular velocity uj ~ , the effects of some forces orig¬ 

inating from it will average down over many rotations; 
the components of the force lying in the plane of rota¬ 
tion are much less effective at producing a net displace¬ 
ment. For uiT ^ 1, the displacement of the PM will 
be aT/oj. Finally, the PM may be subject to stochastic 
forces. If the timescale of the force’s fluctuations is much 
shorter than the coast time, the RMS displacement is 
~ where a = <JaV^, Ca is the standard deviation 

of the stochastic acceleration, averaged over a time pe¬ 
riod T long enough that successive averages are indepen¬ 
dent while still shorter than the coast time. The types of 
forces, their resulting displacements, and the constraints 
required to keep their displacements below 1 m are shown 
in Table m As we will discuss below, these constraints 
appear achievable with careful design and engineering of 
the spacecraft; therefore, we will treat the DSN ranging 
as the limiting source of error for the experiment. 

A. Ranging 

The radial position of the spacecraft needs to be peri¬ 
odically measured with considerable precision during the 
entire coast phase. Here we assume that the uncertainty 
of this ranging measurement is dominated by the perfor¬ 
mance of the Deep Space Network (DSN) |3^. The re¬ 
ported range Icr accuracy is 1 m. The DSN can also per¬ 
form a Doppler measurement of the spacecraft’s velocity 
with an accuracy of 0.1 mm/s; however, successive rang¬ 
ing measurements determine the average velocity more 
precisely, so Doppler information would contribute lit¬ 
tle additional information to the fit. For this reason, we 
only consider ranging measurements in the rest of this 
paper. Negligible uncertainties of 10“^® yr/yr are to be 
expected in the time-base for such measurements, using 
conventional atomic clocks. We note here that the tra¬ 
jectory of the coast, discussed in Sec. m affects the 
telemetry and ranging. While a single DSN antenna can 
be used for the polar coast considered here, a trajectory 
in the ecliptic plane would require a network of ground 
based antennas during the later stages of the coast once 
the round-trip travel time of the DSN signal exceeds a 
few hours. 

The two-spacecraft formation slightly complicates the 
range measurements. However, the short distance be¬ 
tween the DF and the communication spacecraft should 
be easily measurable with negligible uncertainty com¬ 
pared to the Earth-spacecraft distance. An interval of 
~ 3 weeks between measurements, during the entire 
coast, is sufficient to achieve the accuracy required. 


B. Thermal forces 

Temperature differences on the inner surface of the DF 
cavity lead to a net force on the PM. Fairly large ATs 
should be expected from internal components of the DF 
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Source 

Random? Type Magnitude 6R 

Design Constraint 

Ranging 

yes pos. 1 m 1 m 

— 

Thermal gradient 

—External 

—Internal 

no R~^ accel. oq = 

no rot. accel. a = at 

me oj 

AT < 3 X 10-5 K(50o 

AT < 50 K(500 

Charging 

—Electrostatic 

—Lorentz 

(y:q,d accel. a a; . 

(xq,B accel. a = ^at^ 

9max ^ ^ lO^e 

g'max < 3 X lO'^e(^-g^) 

Self-gravity 
— Dipole {x, y) 

—Dipole (z) 

— Quadrupole (Q) 

no rot. accel. a ~ ^ 

no accel. a a; |t^(5(asin0) 

cx d accel. a ~ f^Md ^ i ^2 

r% rc 2 

AM < 0.5 kg 
S{AMs\nQ) < 6 X 10-® g 
/ < 1 X 10-“ 

Residual gas 

yes accel. d 

p < 30 Pa(“° amu)l/2(300 K)l/4 

Solar System uncertainties 

no accel. see text 

— 


TABLE II. Statistical and systematic errors constraining the sensitivity of the experiment. “Random?” is either ‘yes’ for purely 
statistical errors, ‘no’ for purely constant systematics, or a list of potentially random quantities that the error depends on. (We 
conservatively assume that these values are constant and set to their maximum value.) “Type” is either ‘pos.’ for an error 
in the position measurement or ‘accel.’ for a force acting on the PM, optionally with ‘rot.’ to indicate that force is rotating 
(and hence averaged down over many rotations of the shell). “Magnitude” is the parametric size of the effect, and 5R is the 
displacement of the PM due to the effect after a time t. “Design Constraint” gives the required size of various parameters in 
order to keep SR < 1 m after the entire 7 yr coast, assuming the gross parameters given in Table |I] Each source of error is 
discussed further in the text. 


Type 

Displacement 

Constraint 

constant 

\aT^ 

a < 4 X 10-“^ m/s^ 

1 

to 

aoRoT 

V 

oo < 2 X 10-“5 m/s^ 

rotating 

aT 

(jj 

a < 3 X 10-® m/s^ 

stochastic 

~ dT3/2 

a<3x 10-13 “ 

s-^vHz 


TABLE III. Summary of the types of forces acting on the 
proof mass, along with the resulting displacements and the 
constraint on their sizes obtained by requiring that the dis¬ 
placement remain below 1 m for the entire coast. 


spacecraft, including electronics and power sources (espe¬ 
cially if powered by an RTG). However, external sources 
of temperature differences are more constraining because 
they do not rotate with the shell; internal sources are 
discussed at the end of this section. 

The flux of solar radiation incident on the shell of the 
spacecraft leads to a thermal gradient across the outer 
surface of the shell. This gradient can propagate to the 
inner surface of the shell, leading to a net force on the 
PM proportional to (T -|- AT)^ _ ^ 4T^AT, where T 

is the average temperature of the surface and AT is the 
temperature differential between the hot and cold sides. 
Although the shell is rotating, the equilibrium tempera¬ 
ture difference and hence direction of the net force does 
not corotate with the shell but rather maintains its align¬ 
ment relative to the Sun. This leads to a constraint on 
the (non-corotating) inner surface temperature gradient 


of AT(1 AU) < 3 X 10-5 k( 22^)3, using the fact that 
the gradient due to the solar flux falls with distance as 
R~^. Assuming the total mass, size, and rotation values 
from Table |T] and using the heat capacity of mm-thick 
aluminum as a typical value, the temperature difference 
on the outer surface of the spacecraft at i? = 1 AU will 
be on the order of 1 K. Therefore the outer and inner sur¬ 
faces of the shell must be thermally decoupled by some 
insulating layer in order to prevent a large thermal force 
on the PM. 

Such a small inner AT can in fact be achieved 
by adding an insulating layer (e.g. vacuum with low- 
emissivity coatings) between the outer and inner surfaces 
of the shell (see Fig. [^. The temperature difference on 
the inner surface is then suppressed relative to that on 
the outer surface if the shell rotates with a period much 
shorter than the characteristic time of the heat transport 
between the surfaces. (Note that the inner and outer 
shells can be mechanically fixed together and corotat¬ 
ing; it is the rotation of the shell relative to the Sun 
that smooths out the temperature.) The heat transport 
equations are analyzed more fully in Appendix The 
main result is that in the limit of fast rotations, the 
suppression from a single layer of vacuum insulation is 
ATinner = ) where uj ~ 0.6/s is the angular fre¬ 

quency of the shell’s rotation and Trad — 400 s x (^^) for 
1 mm-thick aluminum at 300 K; e is the emissivity of the 
surfaces. A large enough suppression may be generated 
either by making e small or adding a second insulating 
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FIG. 3. Cross-sectional diagram of the rotating shell with a 
layer of vacuum insulation. The surfaces facing the vacuum 
are assumed to have emissivity e. While the material of the 
shell rotates at uj, the direction of the inner and outer tem¬ 
perature differences AT are fixed relative to the Sun. 


layer, as the suppressions from additional layers multiply. 

Temperature differences much larger than that created 
by the Sun will arise from sources internal to the shell, 
such as electronics and power source; unlike those from 
the Sun, these differences rotate with the shell and are 
therefore much less constrained. The temperature dif¬ 
ferences in the plane of rotation can be as large as 50 K 
without significantly displacing the PM. The effect of the 
temperature difference along the axis of rotation is sup¬ 
pressed by choosing the axis of rotation perpendicular to 
the Sun-spacecraft axis, but in the end will be absorbed 
by fitting to the transverse acceleration qt of the PM, 
which is discussed further in Sec. iniDi 


C. Charging 

Electromagnetic forces also affect the PM as it charges 
from cosmic ray collisions. A charged mass inside a spher¬ 
ical (or cylindrical) conducting cavity experiences an ac¬ 
celeration if the mass is offset from the center. Addition¬ 
ally, without magnetic shielding, the PM will experience 
a Lorentz force from (at minimum) the ^nT magnetic 
field in the Solar System, as well as any stray fields from 
the electronics in the shell. If the Solar System magnetic 
field is dominant, the magnetic and electrostatic effects 
have a similar strength with the electrostatic effect giv¬ 
ing a slightly stronger constraint of (/max 1 x lO^e. As¬ 
suming a charging rate of 10 protons/s (similar to that 
expected for the PM in EISA’s drag-free system [H]), 
discharging the PM using UV LEDs [351133] once every 
three days is sufficient to keep the resulting electromag¬ 
netic forces sufficiently small. Most likely, by discharging 
more frequently, the stray fields from the electronics will 
not need to be engineered or shielded all the way down 


to ~nT, but a detailed analysis of the electronics and 
associated fields is beyond the scope of this study. 


D. Self-gravity 

Gravitational interactions between the shell and PM 
are not corrected by the drag-free setup and are there¬ 
fore a potential source of error. To gain an understanding 
of the typical size and type of gravitational force that can 
arise, we consider a generic mass distribution p(r) for the 
shell and assume that p{r) = 0 for |r| < rc, the inner ra¬ 
dius of the cavity in which the DE system operates. The 
PM is taken to lie at r = 0. Then we can perform an in¬ 
ternal multipole expansion of the gravitational potential, 
so that 


cti = Pi + Qijrj + higher multipoles, 

P,=G I dPr' 


r,/3 ’ 


Q^J = \g I Pt 






( 2 ) 

(3) 

(4) 


where we have included only the dipole Pi and 
quadrupole Qij because the higher order terms are neg¬ 
ligible: the maximum size of the acceleration from i-th 
multipole moment is , where d is the PMs 

displacement from the center of the cavity. For the pa¬ 
rameters assumed in Table jlj the £ = 3 and higher terms 
are naturally small enough to ignore. 

Since the shell is rotating, the multipole moments vary 
in time. Expressed (as we have) as tensors, they trans¬ 
form like tensors under a time-dependent rotation ma¬ 
trix. However, the PM displacement does not rotate, 
so the force on the PM from Eq. ([^ has a complicated 
(non-tensorial) transformation. The largest effects come 
from time-independent forces, so we first consider those. 
Each multipole moment has 2.^ -|- 1 independent compo¬ 
nents and one of those (which can be identified as the 
m = 0 part) does not transform under rotations in the 
x-y plane. For the dipole this is simply P^, which leads 
to a constant force perpendicular to the plane of rotation. 
Since we have chosen the plane of rotation to contain the 
Sun, this force does not have a leading-order effect on the 
PM’s motion. However, the plane of rotation must be 
carefully chosen to ensure that the effect remains small, 
and the second-order effect from misalignment is still sig¬ 
nificant. To address this, we leave the constant transverse 
acceleration ot — JA as a free parameter in the ht; then, 
as long as the spacecraft is pointed accurately, the ar 
measured from the fit can be combined with the point¬ 
ing of the spacecraft to accurately subtract this system¬ 
atic. Fitting to Qt accounts for all constant, non-rotating 
forces along the axis of rotation, including the internal 
temperature difference discussed in Sec. |HIB| The coor¬ 
dinate system and the geometry of the situation is shown 
in Fig. [^ 











E. Residual gas 



FIG. 4. Coordinate system used in Sec. |mD] for the rotat¬ 
ing shell. 6 is the misalignment angle between the plane of 
rotation and the direction to the Sun. 


For the quadrupole, the time-independent compo¬ 
nent Q appears on the diagonal: Qxx,yy = Q + 
time-dependent terms, Qzz = —2(3. The magnitude of 
the resulting acceleration (in the x-y plane) is a = Qd 
where d is the displacement from the origin. The bound 

on the size oi Q < is much larger than can be ig- 
’'c 

nored, but is also very conservative, so a more careful 
argument is required. In particular, the bound assumes 
all of the shell’s mass is concentrated near the DF cavity, 
while a more realistic mass distribution would have most 
of the mass located towards the outer surface of the craft, 
changing the parametric dependence of Q. For a uniform 
mass density p, we have Q ^ Gp In ^ ^ In ^, giv¬ 
ing a ~ 3 X 10“^^ m/s^ for the parameters in Table 
Precise engineering of the spacecraft shell, on the order 
of 10“^ fractional precision on the mass and position of 
components, is required to reduce this acceleration to a 
negligible level. The actual form of the spacecraft shell 
could be based roughly on a zero-Q geometry (such as 
a uniform-density sphere or cylinder with equal length 
and diameter) as a starting point and then be shimmed 
to reduce Q to the necessary level. Furthermore, this 
estimate of the required precision is still conservative 
because, with periodic corrections to the offset d pro¬ 
vided by the DF system, the displacement will grow more 
slowly than as assumed here. This argument is only 
intended to estimate the required precision and establish 
plausibility; a detailed analysis of the mass distribution 
and associated uncertainties will be required for any pro¬ 
posed implementation. 


The relevant time-dependent components are Px and 
Py, but these components are naturally small enough 
given the stringent constraint on Q. 


The spacecraft shell will be under a high but finite vac¬ 
uum. Residual gas will induce a stochastic acceleration 
of the PM via molecular collisions. We assume the gas to 
be in thermal equilibrium with the shell and PM, with 
a temperature T, pressure p, and typical particle mass 
p. Since the size of fluctuations of the average acceler¬ 
ation depends on the averaging timescale, the relevant 
quantity to calculate is a = aa\/T where d is the average 
acceleration over timescale t and tra denotes the standard 
deviation of d. Then, for r long enough that successive 
as are statistically independent, the RMS displacement 
of the PM will be given by ~ aT^/^. For the case of gas 
in the molecular regime, we can take r to be the mean 
time between collisions and Ua to be typical acceleration 
imparted by a single collision over time r. Then it is clear 
that (Ta — ^ and r ~ {nAv)~^. For an ideal gas n = ^ 

and so altogether d ~ {pkTY/. The 

constraint on the size of a stochastic force from Table IIIII 
translates into a limit on the pressure p < 10 Pa for 10 
amu gas particles at 300 K. 


F. Solar System objects 


Objects in our Solar System exert irreducible gravita¬ 
tional forces on the PM. Hence, they must be included in 
the model of the spacecraft’s trajectory. Most objects, in¬ 
cluding all of the inner planets, are well-measured enough 
to include in the model with no effect on the experiment’s 
sensitivity [231 HI] • However, the effects of the Sun and 
the Kuiper Belt are uncertain enough that they must 
be included as free parameters in the fit. (As a side ef¬ 
fect, this experiment will provide the best measurement 
of the Kuiper Belt’s mass distribution; see Sec. IV Dj ) 
The Sun is always included as a free parameter since for 
large Yukawa scales A the dominant effect of the new 
force is to unobservably rescale the mass of the Sun. 

Of particular importance is the modeling of the Kuiper 
Belt. Unlike a point-mass-like planet, the mass distribu¬ 
tion of the Kuiper Belt is highly uncertain. In particular, 
a trajectory which passes directly through the Kuiper 
Belt in the ecliptic plane is subject to a large system¬ 
atic uncertainty which is difficult to characterize. One 
could still set a conservative limit using the data from 
such a trajectory, but the resulting limit would be much 
worse than existing limits from planetary motion. How¬ 
ever, a polar trajectory does not pass through the Kuiper 
Belt and allows for a systematic expansion of its gravi¬ 
tational potential (see Appendix . As a result we can 
include a finite number of parameters in our Solar Sys¬ 
tem model. Under the assumption of a very nearly polar 
orbit (so that the Kuiper Belt can effectively be averaged 
azimuthally), we find that the necessary parameters to fit 
to are (equivalent to) the mass, radius, and offset from 
the ecliptic plane. 
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G. Statistical method 


To translate the periodic range measurements obtained 
from the DSN into a measurement of (or limit on) new 
Yukawa-type forces, we must fit the measurements to a 
model that includes both the hypothetical new force and 
all of the known systematic effects. We will consider 
the simplest case where all of the ranging measurements 
Xi are independent and Gaussian with equal variances. 
This case is realized if all of the design constraints in Ta¬ 
ble |TT] are satisfied, although random forces can also be 
accounted for by the fit procedure with a slight gener¬ 
alization. In this case the best-fit (maximum-likelihood) 
parameters of the model are determined by minimizing 






(5) 


Here 9 = {9a} is a vector of the free parameters of 
the model ^{t;9), which gives the expected range of the 
spacecraft as a function of time and the model parame¬ 
ters. The standard deviation of each range measurement 
is (T. We indicate the best-fit parameters by 9. Then 
the uncertainty of each of the 9 is encoded in the matrix 
inverse of the second derivative of 


{v-\t = 


1 9V(g) 

2 d9ad9}, 


( 6 ) 


where the inverse on the LHS is a matrix inverse. Vat is 
the covariance matrix of the best-fit parameters, so that 
the one-sigma uncertainty in 9a is \JVaa- Of course, this 
value depends on all of the derivatives of x^ via the ma¬ 
trix inverse, and the best-fit parameters may have sub¬ 
stantial covariance. 

The expected sensitivity of the experiment can be com¬ 
puted approximately by taking the expectation value of 
Eq. (|^, leading to 


i 


_2 dn{ti) dfijti) 
d9a d9b 


(7) 


also known as the Fisher information matrix. Note that 
taking the expectation value of Eq. is equivalent to 
neglecting terms oc {X^ — , which are suppressed 

in our case since the A6* corresponding to A/i ~ a are well 
within the linear regime, and so inverting Eq. Q is an 
excellent approximation to the true expected sensitivity 
(1A&). Furthermore, computing Eq. 0 involves only first 
derivatives of /x, which is a substantial simplification since 
/i must be numerically integrated with very high precision 
for each value of 9. 

It must be noted that the confidence interval obtained 
from this procedure is only as trustworthy as the model 
itself, so care must be taken to include all relevant con¬ 
tributions. In particular, a small value of x^ at the min¬ 
imum is not a sufficient condition for the validity of the 


A (m) 



FIG. 5. Effect of fit parameters on experiment’s sensitivity 
to a new Yukawa force. The lowest curve results from fixing 
all parameters except GMsun and q; from bottom to top we 
add as free parameters the initial position and velocity and 
transverse acceleration, mass of the Kuiper Belt, and shape 
of the Kuiper Belt as described by its radius and offset from 
the ecliptic plane. 


confidence intervals—one-sigma fluctuations of x^ are of 
size Vs A ~ 14 for N ~ 100 measurements, comparable 
to a l^-sigma systematic effect. Instead we use Ay^ < 1 
as the criterion for exclusion from the fit. 

The free parameters in our simplified Solar System are 
the masses of the Sun and the Kuiper Belt; the Kuiper 
Belt s and ^ moments (defined in Appendix equiva¬ 
lent to the radius and offset from the ecliptic plane); the 
initial position and velocity of the spacecraft; and the 
transverse acceleration of the spacecraft along its rota¬ 
tion axis. The model is then augmented with the param¬ 
eters of the new physics model under consideration, for 
instance a for a new Yukawa force with some scale A, or 
m~^ for a theory of massive gravity. 

The effects of the various fit parameters on the exper¬ 
iment’s sensitivity are shown in Fig. [^for the case of a 
new Yukawa force; the results are discussed further in 

Sec.IVSl 


IV. SENSITIVITY 

A. New forces 

To compute the sensitivity of the experiment to new 
Yukawa forces, we compute the expected inverse covari¬ 
ance matrix according to Eq. 0. modifying the gravita¬ 
tional potential in accordance with Eq. 0 and including 
the free parameters of our simplified Solar System as well 
as the strength a of the Yukawa correction. We present 
the expected 2cr limit on |a| as a function of the scale A 
in Fig. 

At scales A much greater than the distance scale of any 
particular experiment, the Yukawa correction to Eq. Q 
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A (m) 



A(AU) 

FIG. 6 . Expected 2a exclusion limit on the strength a of a 
new Yukawa-type force with range A. The blue solid curve is 
for conservative 1 m ranging precision, while the red dashed 
curve is an optimistic experiment with 10 cm ranging as the 
limiting uncertainty. The shaded region is excluded by plan¬ 
etary tests of Kepler’s third law. 


becomes approximately ^ (1 ~ x + ^ + 

C1(A“^)). The first term is degenerate with rescaling the 
source mass, and the second term is an unobservable con¬ 
stant shift in the potential, so that the leading observable 
correction is a constant radial acceleration a = . As 

a result the limit on |q!| is generically proportional to A^ 
at large A. At scales A > 100 AU our expected limit, as¬ 
suming 1 m ranging, is uniformly two orders of magnitude 
stronger than the best existing limit (planetary tests of 
Kepler’s third law). Ten-centimeter ranging, which may 
be possible in the future [15] , would improve the limit by 
another order of magnitude if it remained the limiting 
source of uncertainty. 


B. Long-distance modifications of gravity 


While Yukawa forces are a generic possibility for new 
physics at long characteristic distance scales, it is gener¬ 
ically difficult to construct complete models of Yukawa 
forces with scales longer than the Earth-Moon distance 
that our proposed experiment is sensitive to without run¬ 
ning afoul of limits on equivalence principle (EP) viola¬ 
tion. One class of models (including DGP gravity |7| and 
some theories of massive gravity m) that respects EP 
while also modifying gravity at long distances relies on 
the Vainshtein mechanism m to screen the scalar mode 
of the graviton at short distances, resulting in a non¬ 
linear power law correction to the effective Newtonian 
potential |i8H5T] . To accommodate a variety of mod¬ 
els simultaneously, we parametrize the correction to the 
gravitational potential as 




GM 

~r' 



( 8 ) 



FIG. 7. Expected 2a lower limit on the graviton Compton 
wavelength m~^ of a new Vainshtein-type power-law contri¬ 
bution to the gravitational potential. For a particular theory 
with index v, the experiment will be sensitive to the region 
below the curve. The vertical lines at | and v = 2 cor¬ 
respond to DGP gravity and ghost-free A 3 massive gravity 
respectively, and the horizontal line indicates the present-day 
Hubble scale which is the scale of interest for models that 
attempt to address the cosmological constant problem. The 
solid blue curve is for 1 m ranging precision, while the dashed 
red curve is an optimistic scenario with 10 cm ranging as the 
limiting uncertainty. The shaded area is excluded by lunar 
laser ranging. 


where v is the “index” of the power law (54/ oc R'') and 
i?* is the Vainshtein radius, the distance below which 
the scalar graviton is screened. Equation Q holds in the 
screened regime, i? <C i?*. In terms of the graviton mass 
m, i?* = (In this section we set c = = I 

so that represents a length scale.) 

For DGP gravity, v = ^ and m~^ ~ is the 
DGP length scale expressed in terms of the 4D and 
5D Planck masses, Mpi and M 5 ; Equation Q reduces 
to 541 ~ rri'jGMR, in agreement with Ref. US). A 
different power law potential with u = 2 results in 
the Vainshtein-screened regime of “ghost-free” (or “A 3 ”) 
massive gravity da SSI; in this case Eq. Q reduces to 
54' ~ m^R?. In both DGP and massive gravity, the value 
of the graviton Compton wavelength relevant to the 
cosmological constant problem is today’s Hubble scale, 
0(iF,("^) ~ 4 Gpc, corresponding to a graviton mass of 
m ~ 10“^^ eV. 

The expected sensitivity of our proposed experiment to 
m~^ as a function of u is presented in Fig. The exper¬ 
iment is sensitive to length scales below the curves; the 
index v characterizes a particular theory. (Here we com¬ 
pare to lunar laser ranging since it is a direct probe of the 
potential in the screened regime; Ref. [S^j has considered 
indirect effects of a graviton mass on weak leasing. They 
do not directly constrain the parameter space in Fig. 
because they rely on the unscreened regime of massive 
gravity where the potential cuts off exponentially.) For 
the particular case of DGP gravity ( 1 / = |), our expected 
2(7 limit on m~^ is 1 Gpc (10 Gpc) for 1 m (10 cm) rang- 
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ing, compared to the current best limit of ~ 100 Mpc 
from lunar laser ranging [HIIllEa]. Thus, our proposed 
experiment is sensitive to the most interesting range of 
DGP parameter space. For steeper power laws (such as 
the ly = 2 potential of A 3 massive gravity; see Fig.[^ our 
experiment becomes even more sensitive relative to lunar 
laser ranging due to the longer baseline of our measure¬ 
ments. 


C. Modified Newtonian dynamics 

Modified Newtonian dynamics (MOND) invokes de¬ 
viations from Newton’s second law at small accelera¬ 
tions a oq — 10 “^*^ m/s^ as an alternative explanation 
of galactic rotation curves that does not require the ex¬ 
istence of dark matter [TBl [19]. However, in this form, 
MOND has a free functional parameter that interpolates 
between F = ma for a ^ oq and F oc for a ^ oq. The 
smallest acceleration (due to the Sun) that the space¬ 
craft experiences during its coast is a ~ 10 “® m/s^, so 
we are always in the deeply Newtonian regime of MOND, 
where predictions are highly dependent on the interpo¬ 
lating function. 

However, there are various indirect effects of MOND 
(due to the external field of the galaxy or the as- 
phericity of the Solar System [55]) that are relatively 
model-independent. These should manifest themselves 
in the Newtonian regime as an extra quadrupolar poten¬ 
tial (5'k ~ Q 2 R^■ Interpreting the limit on a = 2 power 
law from Fig. as a quadrupole contribution leads to 
an expected limit of Q 2 10 “^^ s“^, roughly an or¬ 
der of magnitude smaller than predicted by Ref. for 
most interpolating functions. (It should be noted that 
the precise limit on Q 2 depends on the orientation of the 
spacecraft trajectory, e.g. relative to the Galactic center 
for the external field effect.) Due to these effects, our 
experiment would be sensitive to (or able to constrain) 
broad categories of MOND models. 


D. Kuiper Belt measurements 

As the precision of our proposed experiment requires 
fitting to several parameters describing the Kuiper Belt, 
the experiment also provides a measurement of those 
parameters. Computing the expected sensitivity (via 
Eq. 0) for 1 m ranging in the absence of any new physics 
contribution results in a precise measurement of GMkb 
with an absolute uncertainty of ~ 5 x 10“^GMEarth, cor¬ 
responding to a 0.5% precision if the measured mass is at 
its upper bound of roughly O.lGMEarth [iS]- Addition¬ 
ally, the mass-weighted mean radius and offset from the 
ecliptic plane could each be measured with a precision 

of about 0.2 AU ^ 0 igme^ th )' These would constitute 
the first direct measurements of the mass distribution of 
the Kuiper Belt. Combined with optical and infrared 


observations [SHIM], knowledge of the mass distribution 
could constrain the number and distribution of the small¬ 
est Kuiper Belt objects. 

V. CONCLUSIONS 

We examine the possibility of a space mission to 
100 AU dedicated to the precision study of the 1/i?^ be¬ 
havior of gravity. Such an experiment would extend the 
long-distance edge of our knowledge of the gravitational 
force. 

The ~100 AU baseline enables a much more powerful 
probe of long-distance modifications of gravity because 
their effects, relative to ordinary 1 /i?^ gravity, increase 
with distance. Even assuming a simple new Yukawa-type 
force, the sensitivity at the longest distance scales would 
improve by two orders of magnitude over current limits. 
But in fact, current Yukawa limits are an extrapolation 
from shorter-distance tests and the experiment described 
here would be the first direct test of gravity that actually 
reaches ~100 AU. This is particularly important for val¬ 
idation of various theories of modified or massive gravity 
that provide alternatives to dark matter or a cosmological 
constant; in particular, our proposed experiment is sensi¬ 
tive to the interesting parameter space of the DGP model 
that modifies gravity at the Hubble (~Gpc) scale and to 
the effect of the Galactic gravitational field in MOND. A 
space mission of this type and range is the longest-scale 
direct test of gravity achievable in the foreseeable future. 
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Appendix A: Kuiper Belt parametrization 

The Kuiper Belt (KB) is a collection of small ob¬ 
jects orbiting the Sun beyond Neptune. While optical 
and infrared surveys constrain the luminosity distribu¬ 
tion of the KB, its overall mass is dominated by the least- 
luminous objects which are not detected by these surveys 
and is quite poorly constrained [STIISHj. Despite its small 
mass (less than about O.lMEarth — 3 x 10“^Msun [56| 1 
the large uncertainty in its mass and spatial distribution 
greatly limits our sensitivity to long-range effects. 

To parametrize the gravitational effect of the KB on 
the PM, we perform a systematic expansion of its gravi- 
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tational potential assuming that the spacecraft does not 
pass through the KB itself (i.e. assuming a trajectory 
somewhat above the ecliptic). We will find later that 
a very nearly polar trajectory is required for a reason¬ 
ably simple model of the KB. To perform the expansion, 
change to cylindrical coordinates (s, (p, z) with a fixed ra¬ 
dial offset a: 


elude only Mqo, Miq, and Mqi as fit parameters. (Float¬ 
ing these parameters is equivalent to floating the total 
mass, mean radius, and z-offset of a thin ring represent¬ 
ing the KB.) 

Appendix B: Thermal gradients on a rotating shell 


X = (a -I- s) cos (p 

y = {a +s)sm(p . (Al) 

z = z 

In these coordinates the gravitational potential of the KB 
is given by 


The shell of the spacecraft receives a large flux of 
energy from the Sun, creating a temperature profile 
Touti0,t) on the outer surface which is propagated to 
the inner surface via radiative transfer through a layer 
of vacuum insulation. The full heat transport equation, 
including conduction around the inner layer, is 


'I'KB(r) = —G 


dcp’Az’As’ 


{a + s')pkb{sPzPP’) 
|r - r'l 


(A2) 


where pkb is the mass density of the KB. In analogy with 
a multipole expansion, we now assume s',z' <C |r — r'| 
over the entire support of pkb • In principle, one can also 
expand pkb into a Fourier series in cp to obtain a com- 
plet^series expansion for 4 'kb 7 with coefficients indexed 
by the number of powers of s', z', and appearing in 
the associated term of the integral. 


However, we can greatly reduce the number of rele¬ 
vant parameters with two considerations. First, the radar 
ranging measurement is most sensitive to forces aligned 
with the Earth-spacecraft direction; for near-polar trajec¬ 
tories, the effects of displacements in the ecliptic plane 
are suppressed by ^ which is only a few percent or 
less by the time substantial displacements from the KB 
are accumulated. This fact allows us to consider only the 
radial force, and hence the potential only along the tra¬ 
jectory itself. Secondly, if the trajectory is exactly polar, 
the non-constant Fourier modes of pkb do not contribute 
to the potential: 


^'kb(?'z) = -G 


, +s') Jdp' pkb{s',z',P') 

" " ^{a + s')^ + {r-z')^ ’ 

(A3) 


so only the azimuthally averaged pkb contributes. The 
first few terms are 


^'kb(?'z) = -G 


Moo 


r{rMio + aMoi) 
a(a2-hr2)3/2 + 



where = J dpdzds as'z^ pkb(s, z, p). The effects of 
the i + j > 2 terms are small enough to ignore, so we in- 


T(0, t) = r,ad(T;ut(0, t) - r(0, t)) + Feond ^"^^7 

(BI) 

where F^ad = (^^llfcpd ^’^d Fcond = ^ are the radiative 
and conductive transfer rates. We have assumed that 
the temperature fluctuations are small compared to the 
average temperature T so that the radiative term lin¬ 
earizes. Assuming an outer temperature profile of the 
form Tout{0,t) = T -I-AToutfid® (i.e. a temperature 
gradient that rotates around the shell at angular veloc¬ 
ity w), and an inner temperature profile of the form 
T{9,t) = f + o(t)ed®“‘^*) (so that a{t) is the complex 
amplitude of temperature oscillations on the inner sur¬ 
face), we obtain 

ct iujtd — Fj.a(j ( d) F^ond^- (H^) 

In equilibrium, then, a = 0 and we have a = zATout, 
with 


z = pe ^ = 


rad 


Trad Tcond 


(B3) 


leading to T(0) =T + pATout cos{6 + p — ujt) after taking 
the real part. The suppression factor is then 


ATin 

ATout 


brad 

\/ (brad + bcond)^ + ^2 


(B4) 


which reduces to p ~ for fast rotations. 

For multiple layers, one has to solve as many coupled 
transport equations. However, in the limit p <C 1, the 
effect of additional layers is approximately multiplicative; 
in particular, two layers of vacuum insulation give a ~ 
suppression. 


^ Unlike a multipole or Fourier expansion, this basis is actually 
ouercomplete, and hence not orthogonal. This does not make it 


any less useful for our purposes here, but it is necessary to ensure 
that linearly dependent terms are not included as independent 
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